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Abstract 

We deduced the torque applied on a single stretched and twisted DNA by integrating with respect 
to force the change in the molecule's extension as it is coiled. While consistent with previous direct 
measurements of the torque at high forces {F > IpN) this method, which is simple and does not 
require a sophisticated set-up, allows for lower force estimates. We used this approach to deduce 
the effective torsional modulus of DNA, which decreases with force and to estimate the buckling 
torque of DNA as a function of force in various salt conditions. 
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Most polymers are insensitive to torsion because their monomers are linked by single 
covalent bonds around which they are free to rotate. This property is lost when the polymer 
possesses no single covalent bond about which to release the accumulated torsion. Such is 
the case of a DNA molecule with no nicks (no break in one of the strands), thanks to is 
double-helical structure. This particular feature has very important biological implications. 
First, from a structural point of view, twisted DNA provides an efficient way to compact the 
molecule so that it fits into the cell or nucleus. Second, a negatively twisted (underwound) 
DNA may locally denature thus facilitating the accessibility of its bases to a variety of 
proteins (RNA polymerases [ij, regulation factors[2], etc.). On the other hand, positively 
coiled DNA is more stable at high temperature (it denatures less). Thus thermophilic 
bacteria that live close to the boiling point of water have enzymes (reverse gyrases) that 
overwind the molecule. Because the topology of DNA plays such an essential role in the cell 
life. Nature has evolved a family of enzymes, generally known as topoisomerasesj^, Q, ^(the 
just mentioned reverse gyrase Gj is one of them) that control the torsion and entanglement 
of the molecules. Enzymes that translocate DNA (such as FtsK|7],EcoR124I 8], RScjsl) can 
also apply a torque on the molecule as it is moved along. Thus understanding the behaviour 
of DNA under torsion and estimating the torque arising in a twisted molecule has important 
biological implications. 
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Single molecule manipulation experiments [s], Q, 0, Q, [s, 9, 
twist DNA. In these experiments, a DNA molecule is anchored at multiple points (to impede 
its swiveling) to a surface at one end and to a bead used to apply a force and a torque at 
the other. In the case of magnetic traps, a superparamagnetic bead is pulled by the field 
generated by small magnets and twisting is achieved by rotating the magnetic field ll|, Il2 |. 



With this method the angular position of the bead is imposed and one does not control 
the applied torque. This set-up was nonetheless used to measure the twist-stretch coupling 
in a DNA molecule via the rotational drag of a small bead attached to the backbone and 
allowed to swivel to relax the accumulated torsion|l3|. More recently, optical tweezers jl^ 
have been used to apply a constant torque on an anisotropic transparent particle through 
the angular momentum transfer of a polarized laser beam[l5| that also traps the particle. 
The advantage of the optical tweezers set-up is that it allows for a direct measurement of 
the torque applied on the trapped particle (and through it on the DNA). Its drawback is 
that it involves a rather sophisticated set-up which is difficult to use to explore the low force 
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(below IpN) regime that might be more relevant to biological processes. 

In this paper we describe a simple method to estimate the torque applied on DNA by 
measuring the change in extension of a stretched and coiled molecule as a function of force 
and number of turns. The magnetic trap sy stem, briefly sketched above and employed 
extensively in previous experiments 11, 16,[l3, 13], is used to apply a force F on a magnetic 



bead tethered by a single DNA to a surface and also, by rotating the magnets, to twist the 
molecule by a known number of turns, n. The results of such experiments are qualitatively 
easy to understand on the basis of our daily experience with coiling tubes or cords. Consider 
twisting by n turns a rubber tube of torsional modulus C (usually normalized by ksT in 
the DNA context: C = ksTC) held under an applied force F. Initially the torque F will 
increase linearly with n: T = 27TnC/l, leaving the extension I almost unchanged. Past a 
certain number of turns rib, the associated torque Ff, becomes so large that it is energetically 
less costly for the tube to bend rather than to increase its torsional energy: the tube buckles 
and loops to form a 3D structure called plectoneme, that absorbs torsion as writhe. Further 
twisting of the tube, while leaving the torque unchanged, results in formation of ever longer 
plectonemes. Very similar results are observed (see below and FigH]) when coiling a DNA 
molecule: while its extension I varies little for small n it decreases linearly past a certain 



threshold. As recently suggested [19|, from these observations one can deduce the torque 
applied on the molecule. Indeed the free energy JF of a twisted and stretched molecule 
depends on the force F and the rotational angle, 6 = 27in. The mean extension of the 
molecule at a given force is: / = —dJ-'/dF\0 and the mean torque is: F = dJ-'/dO\F- One 
thus readily derives an expression for the mean torque at a given force F(F, 6) from a 
measurement of the decrease in extension with increased coiling dl/d9\p: 

r(F.<,)^W,«)-£dr(§)^_ (1) 

Since the angular rotation is known and the force and change in extension with rotation 
are easily measurable, the determination of the torque difference is reduced to a problem of 
sampling these variables finely enough to estimate the above integral with sufficient precision. 
The integration constant T{Fq,9) is set by the requirement that there is no torque on an 
untwisted molecule: F(F, 0) = (neglecting the small twist-stretch coupling modulus of 

dnaQ, 3)- 

We measured the DNA extension (of total length Iq 5.4/im) in various salt conditions 
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FIG. 1: (a) Variation of the mean relative extension of a DNA molecule 1/Iq as a function of the 
degree of supercoiiing a = n/Lko in lOOmM NaCl. (b) Variation of the torque in DNA as a function 
of a. These curves have been obtained by integrating as explained in the text Eq{T]from points A 
to B and from points C to D in (a) assuming that the torque in the plectonemic regime at a given 

brce (i.e. at points B and C) is constant. The wiggly continuous lines in (b) are results from ref. 
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and for different values of force F and degree of supercoiiing a = n/Lko (the linking number 
Lko ^ 1500 is the number of times the two strands of the molecule wrap around each 
other), see FiglH Such measurements have been described before UJ, 1221]. Briefly at low 
forces {F less than about 0.4pN) the curves are symmetric. The extension is maximal at 



a = and decreases non-linearly for small values of a due to twist fluctuations [23|, |24 |. 
Past the buckling threshold (at a = as), the molecule coils on itself with a constant slope 
(dl/dalp) to form plectonemes or supercoils. As explained above, while below buckling (i.e. 
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when a < as) the torque increases with increased rotation, it is constant in the plectonemic 



regime 



25l | (i.e. when a > ag). For larger forces (F > 0.5 pN) the curve becomes asymmetric 



as for negative supercoihngs DNA denatures before buckhng at a critical torque: ~ 9 



pN nm 



Id, 



261]. For this reason we have computed the torque only for positive degrees of 



supercoiling, although this estimate should also be valid at low forces {F < 0.5 pN) for 
negative supercoihngs. 

To compute the torque at various forces and degrees of supercoiling using EqUl we start 
from the highest stretching force {Fmax = 3.9 pN in the series shown in fig{T](a)) and the 
highest torque state investigated, namely point A at the buckling transition (a = <Js,max) in 
FiglU^a). The value of torque in A (which served as our reference point) is initially unknown 
but will be fixed by the requirement that T{F, 0) = 0. We then compute the torque at point 
D ( F/) for which ctd < o's,max), by first integrating Eq{T] along path AB (a path of constant 
(T = cta but varying force, see Fig{T]^a)) from F^ax to Fb- Taking into account the fact that 
the torques at points C and B are equal Fc = F^), we then calculate the torque difference 
along the path CD from force Fc = Fb to force Fd along a path of constant a = and 
subtract it from the torque difference along path AB. This procedure is of course valid only 
if both points B and C are in the plectonemic regime. This in effect restricts that procedure 
to values of a > 0.02, where we can identify correctly the buckling transition. To evaluate 
the torque at smaller values, we notice that the values of F as a function of a at high forces 
grow linearly with a. We extrapolate the values of F for one of these curves down to a = 
(requiring F(F, 0) = 0). It does not matter which curve is used: the intercept with the 
ordinate at cr = varies by less than 1 pN nm. We then use these extrapolated values to 
infer from the numerical integration of Eq{T] the values of the torque at other forces and 
values of cr < 0.02. The results are shown in Fig{T](b). Various methods to evaluate the 
derivative have been used and found to yield very similar results (a Savitzky-Golay 27] 
five points second order smoothing method was usually preferred as it is less affected by 
noise due to discrete sampling). The results obtained here are similar, see FigIT](b) to the 
results reported on a different DNA molecule using optical tweezers as a means to measure 
the torque (the slightly different values of the buckling torques might be due to differences 
in the DNA sequences or to different ionic conditions). 

From the data in Fig{T](a,b) we can deduce the effective tension and torque on a bare 
plasmid (circular DNA) unwound by 6% (as often found in Naturej2^). The tension is 
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determined by t 



le value of force at which the DNA's extension is zero at o" = 0.06, i.e. 



F ~ 0.42 pN 



29( 1 (the hne passing through the points BC (at F = 0.42 pN) crosses the 



abscissa at a = 0.06). The buckhng torque at this force is obtained from Fig|T](b): F ~ 6 
pN nm = l.bksT. This value is close to the value where DNA denatures at negative 
supercoilings. This may explain the variable sensitivity of gene expression to sequence and 
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FIG. 2: Variation with force of the effective torsional stiffness of DNA in lOOmM NaCl. The 
continuous curve is a best fit of the high force data to the prediction of ref.jSlI, with C = 94nm. 



The slope of the torque vs. supercoiling curves at low values of cr (see Fig{T](b)) yields 
Cs, i.e. the effective torsional stiffness of DNA. As shown in Figl2] decreases with the 
force acting on the molecule which may explain the low values of torsional stiffness initially 



reported in bulk measurements 32] . This was anticipated by Moroz and Nelson 311] who 
ascribed this variation to a renormalisation of the bare torsional stiffness C by torsional 
fluctuations (that become more important at low forces). Their estimate of this effect 
1 — (C / AA) ksT / AF obtained by a perturbation expansion at high forces) is 



(Cs = C 
shown in Fig 2l 
estimates |lO|, |21 



t allows to deduce a value of C = 94 ± 7nm, in agreement with previous 
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While the value of the DNA bare torsional modulus C does not seem to vary much with 
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salt, the buckling torque of DNA Tt, appears to be much affected by the ionic concentration, 
increasing by as much as a factor 2 at low salt concentrations, see Fig|3l The buckling torque 
increases also with the force 29|] with an approximate power law dependence ~ 
Although the precise values of the buckling torque for DNA may depend slightly on sequence 
(for example AT tracks are known to form bends that may buckle more easily) one expects 
the general dependence of with force and salt to be sequence independent for long enough 
DNA's. Clauvelin et al. [34] obtain the buckling torque from the experimental slopes of the 
torque vs. supercoiling curves using an analytically solvable model of plectonemic DNA. The 
values of buckling torque obtained with their method are very close to the values reported 
in Fig. [31 
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FIG. 3: Buckling torque of DNA for different salt conditions (lOmM, 50mM, lOOmM, 200mM, 
SOOmM NaCl, 2mM MgCl2)- Inset: Log-Log plot of the data and power-law fits (continuous 



lines). The average exponent is 0.72ib0.07. Dotted lines: best fits to a recent model of J.Marko 
(see text). 



J.Marko [25!] has recently suggested a heuristic model to describe the behavior of a 
stretched DNA molecule under twist. In his model, DNA molecules in the plectonemic 
regime partition between an unstretched plectonemic supercoil phase with torsional stiff- 
ness P and a stretched and twisted DNA molecule with persistence length A and effective 
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FIG. 4: Slope of the decrease in extension in the plectonemic regime (the slope of the BC segment 
in FiglUa)) as a function of force and in various salt conditions. The continuous lines are a best 
fits to a recently proposed model [25]. 

torsional stiffness Cs- This model is characterized by only three parameters: the DNA 
persistence length A = 50 nm, its bare torsional stiffness C ~ 90 nm and an unknown 
plectonemic torsional stiffness P, estimated to be about 26 nm. The model makes a number 
of predictions on the variation of extension with a and the variation of Th with F that can 
be compared with experiments. While the predictions of Marko's model are in qualitative 
agreement with our observations (see for example the variation of F;, with force in Figj3l), 
the model cannot explain all the data with only three fit parameters (actually only one P, 
since A can be independently obtained from force-extension measurements, and C can be 
deduced as explained from Figl2]). In particular as can be seen in FigHJ the predictions of 
the model for the variation of 81/ 39 as a function of force does not quantitatively fit the 
data. In some sense this is not very surprising since a description of the plectonemic phase 
with a single force independent torsional stiffness P is an oversimplification that does not 
take into account for example the variation of plectonemic radius with force due to entropic 



repulsion 
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